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Abstract 

A model-independed method of studying the forward-backward 
correlations in symmetric high energy processes is developped. The 
method allows a systematic study of properties of various particle 
sources and to uncover asymmetric structures hidden in symmetric 
hadron-hadron and nucleus-nucleus inelastic reactions. 



1 Introduction 

New data from LHC on soft particle production in pp collisions [U EJ [3] 
open a new chapter in the long history of this problem thus reviving some 
questions which were raised already many years ago. One of these questions 
concerns forward-backward correlations which was extensively studied in the 
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framework of specific models [I] and shown to be useful in discriminating 
between various mechanisms of particle production. 

In the present paper we develop a systematic, model-independent method 
to study forward-backward correlations and show that it may be effective in 
addressing the issue (debated since early seventies) of the number and nature 
of quasi-independent sources of particles contributing to particle production 
in various rapidity regions. Specifically, restricting ourselves to symmetric re- 
actions (like pp or Au — Au collisions) we ask what are the contributions from 
symmetric and asymmetric sources and how to measure them experimentally. 
The problem may be interesting since various models differ substantially in 
this respect. 

The simplest hypothesis is to say that there is just one symmetric source. 
This is the case of the Landau hydrodynamic model [5] and its recent mod- 
ifications [6] where particle production is governed by evolution of a fluid. 
Similar conclusion follows from the simple multiperipheral model [7j suggest- 
ing just one symmetric source in the form of the multiperipheral chain. This 
idea was then reformulated in the parton model [8] and in the bremsstrahlung 
mechanism [SJ. One may of course also consider many symmetric sources. 

A more sophisticated possibility, taking into account the colour structure 
of the colliding systems, was formulated in the Dual Parton Model [TO] . 
Here the number of sources depends on energy of the collision and on the 
type of the projectiles. For p — p collisions, at relatively low energies there 
are two asymmetric sources (chains), spanned between a valence diquark 
from one projectile and the valence quark from the other one. For nucleus- 
nucleus collisions the picture is similar, but the number of asymmetric sources 
fluctuates, depending on number of participants in the two colliding nuclei. 
As energy increases, contributions from symmetric chains, spanned between 
the sea quarks and antiquarks, come into play. 

Simpler ideas were also put forward. In the wounded nucleon model [TTJ 
particles are emitted indepedently from the two colliding nucleons thus cre- 
ating two asymmetric sources [12]. A similar idea, applied to the constituent 
quarks and diquarks, was proposed in [T5] . In the Fritjof model [H] there are 
also two sources, essentially two large diffractively produced clusters (each 
one related to one of the colliding hadrons). 

In the present paper we consider only symmetric collisions (e.g. p-p or Au- 
Au). We show that in this case a systematic study of forward-backward cor- 
relations allows to distinguish between the various possibilities listed above 
and to obtain information about some characteristic properties of the sources. 
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Studies of forward-backward correlations in specific models have a long 
history, see e.g. jU [TO]. Our work was mostly triggered by a recent series of 
papers by Bzdak [H]-[T8], suggesting that a strong asymmetric component is 
present not only in p—p [4J and d—Au (12] but also in Au—Au collisions. This 
observation raises interesting questions about the hydrodynamic evolution of 
the quark-gluon plasma believed to be produced at RHIC [T§] . 

In the next section we formulate the problem in terms of generating func- 
tions. In Section 3 the relations which allow to test various hypotheses in 
p — p collisions are given. Symmetric nucleus-nucleus collisions are discussed 
in Section 4. Some comments on recent STAR measurements [20] are given 
in Section 5. Our conclusions are listed in the last section. 

2 General formulation 

If particles emerge from M independent sources, the generating function for 
the particle distribution in a phase-space region G is a product of generating 
functions describing distributions of particles from individual sources. 

Restricting the discussion to rapidity spectra^], we shall consider two 
sources which are asymmetric with respect to to y — y c . m . — and a third 
symmetric one. Their generating functions are denoted by 4>l, 4>r,, 4>c- 

We shall discuss multiplicity distributions in two intervals of rapidity, 
denoted by A L and A R , situated symmetrically with respect to y — 0. 

Consider first the sum [Al + Ar] . The generating function of the multi- 
plicity distribution in [Al + Ar] can be written as 

®(z;w L ,w R ,w c ) = Y. P (n)z n = [<t> L {z)\^[M^T R [M^T c ■ W 

n 

where wl, wr, wc are numbers of the relevant sources. 

Assuming now that the splitting between and A# of particles from 
each source is random (i.e. it follows the binomial distribution) we have for 
the joint distribution in Al and Ar 

<j>m(zL, Zr) = Pm( n L,n R )zl L Z R R = 4> m (pLmZ L +PRmZR) (2) 

n L ,n R 

where pLm and pR m are probabilities that a particle emitted from the source 
1 Our discussion applies to any variable symmteric with respect to an axis or a plane 
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labelled m ends up in A L or in A R . Thus pim +PRm = 1- Symmetry implies 



Pll=Prr=P+; Plr = Prl=P- = 1~P+; Plc = Prc = ( 3 ) 
Consequently 

^K-^L, Zr\ IV L, W R , Wq) = 

= [<t>{p + z L + P -Z R )] WL \<$>{ V _z L + p + z R )] WR [Mzl/2 + z R /2)] wc . (4) 

For symmetric collisions we have wl = w R = w and <Pl{ z ) — <Pr{z) = 
4>(z). We also note that the distribution of particles in one of the considered 
intervals, say Ax,, is evaluated from ^(z) = &(zl = z,z r = 1). 

From these formulae one can evaluate all moments of the joint distribution 
in Al and A R , as well as the moments of the distribution in one of the 
intervals, in terms of the moments of the distributions describing the sources. 

When only symmetric or only asymmetric sources are present, one can 
derive relations between the joint moments (describing the forward-backward 
correlations) in terms of the moments characterizing the distribution in one of 
the intervals. These relations provide demanding tests for these hypotheses, 
allowing to distinguish between various mechanisms of particle production. 
When both symmetric and asymmetric sources contribute the relations allow 
to obtain information on distributions characterizing the sources themselves. 

In the next two sections we discuss some of these relations. 



3 Relations between cumulants 

In this section we derive relations between the cumulants of the joint 
distributions in A + and A_ and the cumulants fi of the distribution in one 
of the intervals. 

In terms of the generating functions §(zl,z r ) and *&(z) the cumulants 
are defined as 

_ d k+l {log<S>(z L ,z R )} _ d*{\og^(z)} 

o z L o z R [ ZZ _ 2JJ=1 ] az [ 2 =i] 

Since logarithm changes the products in &(zl,Zr) and ^(z) into sums, we 
immediately obtain 

fk+i - fkl = \ [p+ +l + P- +l ~ p\p l - ~ P+P-]fk + i = 9k9ie 2 f k+ i (6) 
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where 



2fi 



d^hg^iz)] 2 ™} 



(7) 



are the cumulants of the particle distribution in [A + + A_] coming from the 
two asymmetric sources and 



Note that in ([H]) the dependence on the number of sources and contributions 
from the symmetric sources drop out. 

The cumulants fi and fki can be determined from the standard factorial 
moments F{q = Foi = Fi and from the joint factorial moments in two intervals 



The relevant relations become rather involved at high orders. The first 
few are listed in the Appendix. 

One sees that, generally, the result depends on two functions 0, (fic and 
one parameter (p + or The resulting relations involve not only the ob- 
served moments but also the moments of the distributions characterising the 
asymmetric sources. Thus, if p + 7^ p_ (i.e. if the sources are indeed asym- 
metric), one can - from the measured Fki and Fk+i - obtain information about 
fk + i characterizing the distributions of particles from asymmetric sources. 

If there are only asymmetric sources we have obviously fk+i = fk+u so 
that all quantities entering the relations ([6]) can be measured. In other words, 
in this case (E]) represent identities between the measurable quantities which 
must be satisfied if the symmetric sources are not present. 

The relation for k = I — 1 is of particular interest, as it allows to deter- 
mine the parameter p + = 1 — p_ and thus to determine the distribution in 
rapidity of the particles from a single asymmetric source [15J. This determi- 
nation is of course valid only if the two-sources idea is satisfied by data. The 
other relations can be used to verify this assumption 

When only the symmetric source is present, all moments (pi vanish and 
we obtain f M = f k+h implying 



a really strong constraints. Naturally, identical result is obtained when p + = 
i.e. when the asymmetric sources become symmetric. 







Fu — Fk+i, 



(10) 
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4 Nucleus-nucleus collisions 



Investigation of forward-backward correlations in nucleus-nucleus collisions 
may be of particular interest, as it allows to verify to what extent the asym- 
metric components survive the period of thermalization and hydro dynamical 
expansion which are believed to determine the outcome of the collision. 

To apply our analysis to this case, one has to take into account that the 
numer of sources may fluctuate, depending on the numer of collisions and/or 
wounded nucleons [TTJ [HI (2TJ |22] . Denoting the number of left (right) movers 
by wl(w r ) and the number of symmetric sources by Wc we obtain for the 
generating function of the joint distirbution in Al and A R 

$w(zl,zr)= E W(w l ,w r ,w g )®(z l ,z R ] w l ,w r ,w c ) (11) 

w L ,w R ,wc 

where W(wl, w r , wq) is the probability distribution of the numbers of sources 
and Q(zl, z r ;wl,wr,wc) is given by (j3J. The generating function of the 
distribution in one of the intervals, say A L , is ^ w (z) = § w (zl = z,z R = 1). 

Fluctuations in number of sources imply that the generating, function 
& w (zl, z r ) is no longer a product of functions desribing the sourcesj. There- 
fore relations between cumulants become rather involved. It turns out that 
somewhat simpler relations are obtained for the factorial moments. For sym- 
metric processes, where W(wl, w r , wq) = W(w r ,wl,wc), the simplest ones 
read 

F 2 — Fn = e 2 <L 2 — LiRi>; 
F 3 -F 12 = e 2 {<L 3 -L 2 R 1 > + < d[L 2 - L X R X ] >} ; 
F 4 -F 22 = e 2 {<L 4 - L 2 R 2 >+2 <d[L 3 - L 2 R 4 ]> + <C 2 [R 2 - L 1 R 1 }>}; 
F 4 - F 13 = e 2 {(1 -p+p-) <L 4 > -t 2 <L 3 i? 1 > -3p + p_ <L 2 R 2 >) + 

+ ^e 2 <C x [Lz - L 2 Ri]> +- A e 2 <C 2 [L 2 - L 1 i? 1 ]>(12) 
where < ... > denotes the average over the number of sources and 

_ <p mz)}^} _ i mz)\ wR } n _ {[<t>c{z)r°} 

L " L ~ Tl ' Ui ~ Tl ' 0i ~~ Tl ' { l6 > 

az l [*=i] az l ^=i] az l ^ 

2 Needless to say that, if there are no correlations between numbers of sources, i.e. if 
the probability W(wl,wr,wc) is a product of three factors, the situation reduces to the 
one described in the previous section. 
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are moments of the distribution in + A# produced by left, right and sym- 
metric sources, respectively. They cannot be directly measured, so ( 1T2|) can 
be used to obtain information about them. Such analysis faces, however, a 
difficulty: One sees that the R.H.S. of ffT2j) depends not only on the structure 
of the sources but also on correlation between the numbers of various sources. 
To disentagle these two effects it is necessary to study processes with various 
nuclei and at various centralities |17j . 

If the symmetric sources are absent, (fl2|) simplifies to 

F 2 - F u = e 2 [< L 2 > - < L X R X >]; 
F 3 -F 12 = e 2 [< L 3 >-<L 2 R 1 >]; 
Fa - F 22 = e 2 [< L A > - < L 2 R 2 >]; 
F 4 - F 13 = e 2 [(l - p + p_) < L 4 > -e 2 < L 3 #i > Sp+P- < L 2 R 2 >]. (14) 

Finally, if only symmetric sources are present we return to the simple 
relation (JTUJ) - 

5 Three intervals 

Fluctuations of the number of sources in nucleus-nucleus collisions are diffi- 
cult to control because even precise determination of centrality of the collision 
is not suffcient to guarantee a fixed number of sources [T71 [181 ED [22] • To 
improve this, STAR collaboration measured correlations in three intervals 
[20] . Apart from the A^ and A#, one adds the third interval Ac, located 
centrally around y cm = and not overlapping with A^ and Ar. The corre- 
lations between A^ and Ar are measured under the constraint that a fixed 
number of particles, ric, was found in Ac. Particle multiplicity in Ac is 
obviously related to the number of sources and therefore one may expect it 
to be helpful in estimating this number on event-by-event basis. 

The extensive general discussion of this data was given in |22j . They were 
also analyzed in [17] and [21] within specific models. 

We would like to add three comments. 

(i) The relations derived in previous sections remain intact if one adds the 
condition that a certain number of particles is observed in the central interval 
Ac- This should be clear from the derivation: replacing the probabilities in 
(CQ) and ([2]) by conditional probabilities (fixing the number of particles in Ac) 
does not change at all our argument. 
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(ii) When symmetric and asymmetric sources are present, restricting nc 
has only a limited effect on reduction of fluctuations of the number of sources. 
This can be seen by considering the distribution of particles in Ac. The 
relevant generating function is 

$c(^c) = [0(1 -Pc + PCZ C )] WL+WR [<M1 - qc + qcz c )] wc • (15) 

where pc is the probability that a particle from an asymmetric source lands 
in Ac, and qc is the probability that a particle from a symmetric source 
lands there. 

One sees from (|15[) that the particle distribution in Ac, does not depend 
on the difference W- = wl — wr and thus fluctuations of u>_ are not restricted 
by (|T5|) . Moreover, ([15]) implies that the distribution of ric depends on both 
w + = wl + wr and wc- Since the forward-backward correlations induced by 
fluctuations of asymmetric sources are, generally, weaker than those induced 
by the symmetric ones, its is important to separate the two contributions. 
We conclude that although measurements at a fixed nc may be helpful, 
probably some additional model assumptions are necessary to disentangle 
this problem. 

Let us also note that in [211 122] only the symmetric sources were discussed 
and therefore this aspect of the problem did not appear. 

(iii) Using the methods of sections 3 and 4, relations can also be de- 
rived for the joint moments of the distribution in the three intervals. The 
only difference is that one has to consider the generating function of three 
variables. 

6 Summary and comments 

A systematic, model-indepedent method of studying the forward-backward 
correlations in particle production is developped. It is shown that it pro- 
vides a useful tool for determining the structure of the sources of particles 
created in high energy cllisions. In particular, it may be used to uncover 
left-right asymmetric components present in these processes and to study 
their properties. This point is of interest since, as explained in Introduction, 
various mechanism of particle production differ in their predictions for the 
presence and/or intensity of such asymmetric sources both in {p — p) and 
A — A collisions. 
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The existing data [23] show reasonable agreement with the idea that just 
two asymmetric sources dominate the observed correlations in p — p collisions 
[H HOI [16]. A similar conclusion was obtained recently in the analysis of 
the Au — Au collisions [21], where data could also be explained without 
any symmetric contribution being present [18j . This seems not to be the 
case [22] for the more restrictive STAR data [2ojf|. These conclusions have 
of course important consequences for selecting the possible mechanisms of 
particle production. 

As was indicated in Section 5, the analysis of the heavy ion data may 
require additional information about the correlation between the various par- 
ticle sources. Even in this complicated situation the method proposed here 
can, however, clearly distinguish whether symmetric or asymmetric sources 
dominate the process in question. 

New data from LHC [H [21 [3] show that the multiplicity in the central 
rapidity region increases with energy much faster than expected from simple 
extrapolations of the trends observed at lower energies. One possible expla- 
nation is that, as predicted in some models [TO] , a new symmetric source 
of particles may be excited at these high energies. Studying the forward- 
backward correlations using methods developped in the present paper should 
be helpful in verification of this idea and, possibly, identification of this new 
component, as well as in investigation of its properties. 

A question which may be studied by the methods proposed in this paper, 
is the comparison of the forward-backward correlations in p — p and in heavy 
ion collisions. Such comparative studies at LHC energies would allow to 
obtain information about longitudinal dynamics of quark-gluon plasma, the 
problem which is barely touched by the existing analyses. Forward-backward 
correlations are created at the very early stage of the collision (see, e.g. [22] ) 
and, apparently, survive the evolution of the system. It remains, however, an 
interesting question to what extent they are distored during this evolution. 

Finally, let us emphasize that our approach is not restricted to rapidity 
distributions: it can be used to study correlations in other variables as well. 
One interesting possibility is to repeat the standard analysis in rapidity re- 
stricting, however, the azimuthal angle (keeping of course forward-backward 
symmetry). This may provide useful information on correlations in the di- 
rected flow. 

Another attractive possibility is to consider correlations in tranverse mo- 
3 See, however, [T7|I5I]. 
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mentum at a fixed and/or opposite rapidity. Since the p t distribution in 
minimum bias sample is now becoming accessible in a rather broad range 
Pfl 121 |3], the full potential of the method can be explored. Such investigation 
may help in disentangling the jet structure in the low and medium p t regions. 

7 Appendix 

Some explicit relations between cumulats and factorial moments (for sym- 
metric collisions) are listed below: 

fi = F\\ fn = Fn — i<\ 2 ; / 2 i = F21 — 2Fxi.Fi — F 2 F\ + 2F-J 3 ; 
f 22 = F 22 - 4F 21 Fi - 2Fl - Fl + 4 [F 2 + 2F n ] F? - 6F?; 
hi = F 31 - 3F 2 i/i - F 3 F X - 3F n F 2 + 6[F n + F 2 ]F? - 6F X 4 ; 

f 2 = F 2 -Fl- f 3 = F 3 -3F 2 F 1 + 2Ft, 
U = F 4 - 4F 3 F! - 3F| + 12F 2 F? - 6F?. (16) 
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